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Using preliminary and published Higgs search data we investigate whether any evidence supports 
the possibility that the Higgs may be mixed with other neutral scalars. We combine the positive 
evidence for the Higgs at 125 GeV with search constraints at other masses to explore the viability of 
two simple models. The first Higgs 'friend' model is simply a neutral scalar mixed with the Higgs. 
In the second Higgs 'accomplice' model the new scalar has an enhanced coupling to photons due to 
couplings to additional charged fields. We find that the latter scenario allows improvement in fitting 
the data by accommodating enhanced h — > 77 rates and suppression in other channels for a Higgs 
mass of 125 GeV. Small excesses at other masses allow the additional scalar to further improve the 
fit to the data, particularly if it has mass in the vicinity of 210 GeV. Due to observed event rates 
at 125 GeV and strong limits in high mass Higgs searches, mixing angles 9 > 7r/4 are typically 
disfavored at the 95% confidence level, depending on the mass of the scalar. 

o 

£ I. INTRODUCTION 

I—* 

The recent announcement of significant evidence for the Higgs boson [IJ [2], gathered at the CERN LHC, brings 
particle physics into a new era of discovery. Further data may confirm this signal to be due to the Standard Model 
(SM) Higgs, however there are already emerging, albeit very weak, hints that production cross-sections and/or decay 
rates might not be quite as expected for the SM Higgs. If these hints persist then it may be the case that the discovery 
Qh of the Higgs comes accompanied by convincing evidence for physics beyond the SM (BSM physics). From a theoretical 
perspective some modification of Higgs physics has long been expected, since substantial theoretical motivation for 
BSM physics is aimed at resolving the hierarchy problem, which is concerned with the Higgs sector. Solutions to this 
problem often require the existence of additional electroweak-charged states and/or additional scalars coupled to, or 
l— "™ 1 mixed with, the Higgs. 

I Even if one abandons the hierarchy problem as motivation there is always the possibility that additional hidden 

sectors exist, perhaps related to dark matter. The Higgs sector of the SM contains a super-renormalizable Lorentz 
&\ and gauge invariant operator, which can easily accommodate couplings to new hidden sector physics, the so-called 
'Higgs Portal' [3HT5], Such couplings may allow for Higgs decays to neutral particles, leading to an additional invisible 
width for the Higgs. The main consequence of this scenario is that all detectable branching ratios become equally 
suppressed, leading to a democratic reduction in the Higgs signal rates. In addition to this, hidden sector scalars can 
f — also mix with the Higgs through the Higgs portal interaction. 

Motivated by these simple considerations we study the implications of the LHC Higgs searches on simple models of a 
singlet scalar mixed with the Higgs. Although simplified, we believe these models should map on to some theoretically 
' . I motivated scenarios, and have applicability beyond simply improving fits to Higgs data. 

^ If additional neutral scalars mix with the Higgs the mass eigenstates and interaction eigenstates are not aligned and 

the properties of the Higgs are altered, similar to the so-called 'Higgs look-alike', or 'Higgs Friend' scenario [Pol - EU] . 1 
If this mixing alone is present the individual branching ratios of the Higgs remain the same and the production of the 
Higgs is suppressed. The overall effect is thus, at the level of the current Higgs searches, indistinguishable from the 
case where the Higgs has an additional invisible width. However, with careful study the two can be distinguished. 
For example, an invisible width can be measured by searching for mono-jet signals coming from initial state radiation 
in Higgs production [141 121H25) . On the other hand, mixing with a neutral scalar can be confirmed more directly by 
searching for the extra scalar, exploiting the fact that it inherits many properties of the Higgs. In this work we are 
concerned with the latter scenario. We assume that the Higgs is present, with mass of 125 GeV, and then consider 
limits from the Higgs searches on a Higgs friend. Of course, it is plausible that the friend might be much more massive 
and thus effectively decoupled, in which case it would be very difficult to unambiguously confirm its presence. We 
consider fits to this scenario as well. 
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1 Although we do not assume that the neutral scalars have couplings to extra colored particles, as in |19| . the scenario considered here is 
sufficiently similar in spirit to the proposals in |19| that we adopt the 'Higgs Friend' terminology. 
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As stated, this simple Higgs friend scenario leads to suppression of all Higgs signal rates, independent of the 
particular search channel under consideration. However, early evidence from the LHC suggests that while some small 
degree of suppression in most channels may provide a better fit to the data, the h — > 77 channel, which is driving 
the statistical significance of the discovery, is possibly enhanced to some degree. As such, it is clear that the Higgs 
friend scenario alone will not lead to significant improvement in fitting the data when compared to the SM. For this 
reason we consider the addition of charged vector-like fields, which could be scalar or fermionic, and couple to the 
Higgs friend. 2 This interaction enhances the coupling of the friend to photons at one-loop. Once the Higgs mixes 
with the friend this can enhance the Higgs decays to photons, allowing a better fit to the data than the SM Higgs. 
We call this the 'Higgs accomplice' scenario. 

In order to test the viability of the Higgs friend and Higgs accomplice scenarios it is necessary to confront these 
models with data. Since the first convincing hints of a Higgs were reported in December 2011 there has been significant 
interest in determining how best to extract Higgs couplings from the data, leading to a number of studies [2"oT - (55] . We 
do this by approximating the likelihood functions for the Higgs signal strength in particular decay channels by using 
the best fit values or, in one case, the expected versus observed 95% exclusion contours provided by the collaborations. 
This information can then be used to estimate best fit points and confidence contours for the models considered. 

The finite mass-resolution of the Higgs searches introduces additional subtleties when considering the presence of 
two Higgs-like scalars contributing signal events in the searches. Whenever the mass separation of the Higgs and friend 
is much greater than the mass resolution of the searches we can take the product of the individual likelihoods, since 
the searches at different masses are approximately independent. However, when they are close enough in mass that 
signal events from both scalars cannot be considered separately a more sophisticated likelihood must be constructed. 
Without performing a collider simulation we suggest a crude method by which to construct the combined likelihood, 
which we believe captures the dominant features of the likelihood function for the Higgs and the extra scalar in this 
case. Rather than using this combined likelihood to make precise statements about scenarios where both scalars 
are close in mass, we instead use it to estimate the mass range in which the searches are effectively independent, to 
determine when the combined likelihood can be trusted. We then find the best fit parameters and associated 95% 
confidence contours for both scenarios. 



Before presenting our results in Sec. |TV 
Sec, [il] and our statistical methods in Sec. 
conclusions in Sec. IVII 



we will briefly review the Higgs friend and Higgs accomplice scenarios in 
III We also consider precision electroweak constraints in Sec. [V]and draw 



II. HIGGS FRIENDS AND ACCOMPLICES 



We consider a simple set-up in which an extra field, s, mixes with the neutral Higgs through a Higgs portal coupling. 
In the mass-eigenstate basis the two neutral scalars are h and s, which are related to the interaction eigenstates through 

( h \ _ f cos 9 — sin 9 \ / h \ 

\ s J ~ \ sind cos6> J ' \ § J ' [ ' 

This is the Higgs friend scenario. We also consider the Higgs accomplice scenario in which s couples to additional 
charged particles. At one loop this leads to a coupling 

C = ac hm aF> u 'F ia , , (2) 

where c^ 77 is the usual SM coupling of the Higgs to photons, and a parameterizes deviations from this coupling. 3 
Without loss of generality we impose < 9 < n/2, and if 9 < 7r/4. We will refer to the scalar at 125 GeV, h, as the 
Higgs. 

All relevant Higgs production cross-sections at the LHC now come suppressed by a factor of cos 2 9 and every decay 
width is suppressed by the same factor, with the exception of decays to photons which are now accompanied by a 
factor of (cos 9 + a sin 9) 2 . Since decays to photons are far subdominant then, to a good approximation, all branching 
ratios remain the same as the SM Higgs, with the exception of the branching ratio to photons which is accompanied 
by the factor (1 + atan#) 2 . Thus for the search channels h — > bb, tt, WW, ZZ the total event rate normalized to the 



2 It is often the case that additional electroweak-charged fields are present in extensions of the SM Higgs sector, so the introduction of 
extra charged fields is a plausible augmentation of the Higgs friend scenario. 

3 In general a can either be positive or negative. 
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event rate for a SM Higgs, otherwise known as the strength modifier, fi, is simply 11 = cos 2 9 and for h — > 77 it is 
/i 77 = (cos# + asin0) 2 . 

Production of the Higgs friend, s, is suppressed by a factor of sin 2 9 compared to SM Higgs production. Whenever 
771$ < 2™^, the strength modifier for the friend in the diphoton channel is (sin# — acos9) 2 and is sin 2 9 for all other 
channels. 

Whenever m~ s > 2m^ the trilinear scalar interactions allow for the decays s — > 2h. These decays could lead to 
interesting signatures, such as 46 final states, however such signals are not currently accessible at the LHC. 4 As we 
are considering the sensitivity of the dedicated Higgs searches to Higgs friends and accomplices we can treat this 
additional width as invisible. Given the physical masses, mixing angle, and scalar potential parameters, one can 
determine the magnitude of this interaction, which is essentially a free parameter, and the resultant width (see e.g. 
|13|). As the trilinear coupling is a free parameter we do not lose generality by taking the invisible branching ratio as 
a free parameter. 5 We can express this branching ratio in a model-independent sense as 

1/2 

BR(S -> 2h) = K I 1 - 4— | I , (3) 




where k is the branching ratio in the limit — > 00 and the kinematic factors are included such that the branching 
ratio vanishes at threshold. 

Hence for a Higgs friend or accomplice the strength modifiers become 



^ = sin 2 6» I-kWI-4— I , ^ 77 = (sinA - acos(9) 2 1 - kW 1 - 4— | , (4) 





for the s — > bb, tt, WW, ZZ and s — > 77 channels respectively. 

Counting parameters, we have the recently determined parameter of the SM, mr, along with four new free param- 
eters TOg, a, 9 and k. In light of recent LHC data we set = 125 GeV and consider the mass range 120 < mg < 600 
GeV. As argued above, to a good approximation, the only search channel sensitive to the parameter a is for the 
decays h — > 77. 

One could also consider coupling the scalars to additional colored fields, which would lead to enhanced production 
in the gluon fusion channel. We will not consider this scenario here for two reasons. Since the enhancement of the 
diphoton channel, and suppression of the non-diphoton channels, can be easily accommodated in the Higgs accomplice 
scenario there is little to gain by boosting the Higgs production in this way and the introduction of this additional 
parameter will not lead to a significant improvement in fitting the data. Also, all Higgs production cross-sections in 
the models considered here are re-scaled in the same way, so above 150 GeV, where the diphoton searches are not 
sensitive, one can employ the reported fits for all sub-channels combined. Whereas to consider boosting gluon fusion 
alone means that all production sub-channels should be treated independently and the relative contributions of gluon 
fusion and vector boson fusion should be considered independently. As these relative contributions and subsequent 
likelihoods must be estimated somehow, this leads to the introduction of further error. 6 



III. STATISTICAL METHODOLOGY 



A. Likelihood estimation 



Before considering scenarios with multiple Higgs-like scalars we briefly review some methodology regarding Higgs 
signal strength likelihood functions, which determine the compatibility of Higgs-like signal with the Higgs search 
results. 7 Within a particular search channel, limits on a single Higgs particle of mass rrih are expressed in terms of the 



4 We thank Christoph Englert for discussions on this point. 

5 Both h and s could also have additional widths to invisible states. For s this is automatically accommodated in this analysis since the 
invisible width is a free parameter. For h the overall effect is to democratically reduce event rates. 

6 For a treament of this scenario with the friend decoupled see |35| . 

7 Throughout we employ a frcqucntist approach. Employing a Bayesian approach allows the likelihoods for strength modifiers to be turned 
into probability density functions when one includes priors. In this case, if strength modifiers depend on additional parameters then a 
Jacobian must be used when changing variables, so that the mode of the probability density function may not occur at the same place 
as the maximum of the likelihood (see for example |36|V The non-invariance of Bayesian estimators under reparametrization motivates 
our choice of a frcqucntist approach. 
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strength modifier /i, which relates the signal strength to that of a SM Higgs at a given mass, \i — n s /(n s ) , where 
n s is the number of signal events expected for the particular search channel. Given the number of observed events 
n D b s , one can construct the likelihood function C(n f, s \fJ>, 6) which is a function of the parameters [i and 6. Here 6 
stands for a set of nuisance parameters, which are fitted from the data to account for systematic effects and unknown 
background estimation parameters. 

The standard quantity used to test hypotheses or set limits on /i is the so-called profile likelihood ratio [37] 

L(n obs \jl,0) 

where is the value of that maximizes the likelihood for a specified value of fi, while fi and 6 are the maximum 
likelihood estimators for /i and respectively. We do not have access to the full likelihood functions £i(n f, s |/Lt, Oi) for 
the different channels, but given the information available from the experimental collaborations we can reconstruct 
approximate profile likelihood ratios Ai(/x). In order to combine results from multiple search channels and different 
experiments, one should in principle calculate the profile likelihood ratio from the combined likelihood, which can be 
taken as the product of the different likelihoods if the channels are independent. Given that we do not know the full 
likelihood functions, we will take the combined profile likelihood ratio as the product of the A^(/x) reconstructed from 
single channels or experiments 

i 

It should be noted, however, that this approximation introduces an additional source of error if significant correlations 
between channels arise, possibly through the nuisance parameters. 

In an abuse of terminology, we will henceforth refer to A(/x) as the likelihood for fi. In order to reconstruct Ai(/z), 
we note that, as described in [551 133 13H] > in the limit where the number of events is sufficiently large, with n i, s > 10, 
the likelihood for a given channel can be approximated by 

A( M ) =s exp-^) 2 / 2 ^ , (7) 

where a b s is in general a function of \x. For the 7 TeV run the best fit strength modifier p, and the error cr b s in 
individual 77, tt, bb, WW and ZZ search channels are reported by the ATLAS and CMS collaborations [33 ED] as a 
function of the Higgs mass. For most channels the Gaussian approximation, which assumes that a b s is independent 
of /i, works well and using the best fit parameters and uncertainties the 95% confidence limits can be reproduced well. 
However in the ZZ — > 41 channel the likelihood function is clearly not Gaussian, as can be seen from the asymmetric 
confidence contours in [391 140j . Thus the Gaussian assumption is not valid and its use can introduce artificial bias 
into parameter fits. 

We choose to approximate A(/Lt) as a two-sided Gaussian, since this captures the approximately Gaussian nature of 
the likelihood and employs the three pieces of information available at a given mass, namely the best fit point and 
two values of the log- likelihood away from the best fit point. Using the 7 TeV data we can test this approximation by 
taking a on either side of the best fit value of /i from the la values (CMS [40]), or A(— 2 log A) = 1 contours (ATLAS 
|39| ) provided for the individual channels. Although this approximation is crude, using the 7 TeV data we find 
that when combining all channels and comparing with the reported combined best fit values the two-sided Gaussian 
assumption fares reasonably well, and typically better than the standard Gaussian approximation with symmetrized 
errors. It should be kept in mind that errors of 0(10 — 15%) are typical using this approach, combined with an 
inherent error due to digitization of the data, which we estimate to be as large as 0(10%). 

In some cases /x is the only free parameter, however in more complicated models involving modified Higgs couplings 
or additional invisible decay widths fi becomes a function (in general different for different channels) of the additional 
parameters of the model, where u> denotes all the additional parameters and the superscript denotes the 

particular search channel. One can find best fit parameters by maximizing the likelihood function and, since the 
quantity —2 log A should approximately follow a chi-squared distribution [37], one can also test the hypotheses of 
different models or signal strengths. 



B. Likelihood for multiple scalars 



Thus far we have only been concerned with models in which the hypothesis is of a single Higgs particle. However in 
this work we consider models containing two Higgs-like scalars of mass m~ h and mj, and we must estimate a combined 
likelihood for both. 
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The individual searches have differing mass resolutions, from as small as 1 — 3% in the h — > 77 and h —> ZZ —> 41 
channels up to 20% in the h —> tt and h —> WW channels. Whenever masses are greatly separated, i.e. \m~ h — ms\ ^> <Ji 
the hypothesized signal from one does not contaminate the search for the other, making the searches effectively 
independent. In this case the likelihood can be taken as the product of the two independent likelihoods 8 

Ac [^~ h '^s,^,m fl ,m s ] = A c [/x^,u?,m^] x A c \p, s ,u,m s ,] ■ (8) 

Whenever the Higgs and friend or accomplice are separated by mass splittings within or close to the mass resolution 
of a given search channel, |m^ — mg| < ct;, the situation becomes more complicated. Without performing a full 
simulation we can still estimate the likelihood in such a scenario based on the mass resolutions provided. However, 
in a conservative approach, we will not use this estimate to make precise statements about fits in regions where the 
signal from both scalars overlap, but will instead use it to determine the mass range in which the factorized likelihoods 
for the scalars can be trusted. 

To understand how we estimate the combined likelihood whenever \m~ h — m%\ < <J%, one can first focus on a single 
search channel and consider a hypothetical situation in which signal from a SM Higgs, of mass m^, is present in the 
data. Performing a search for a Higgs of mass nth, with cuts optimized for this mass, one expects to observe a certain 
number of signal events n s (m,h), and to reconstruct a strength modifier at that mass of fi(mh) ~ 1, up to statistical 
and systematic errors. However, due to the finite mass resolution, a certain number of events, originating from the 
Higgs of mass m^, may also pass the cuts for a Higgs search for a different mass m' h . Hence looking at searches for 
different masses one expects to observe a certain number of events n s (m' h ) < n s (m,h) and to reconstruct a strength 
modifier at that mass fj,(m' h ) < /i(m/j), even though the true Higgs mass is rnh- This makes intuitive sense: for a SM 
Higgs at mass m/,, with a finite amount of data one would not expect the reconstructed strength modifier to be a 
precise delta-function but rather it should follow some distribution which is peaked at m^. 

Given that we know /1 cx (n b s — n&) = n s then, regardless of the mass-dependence of the backgrounds, we need 
only know the dependence of the eventual signal on the Higgs mass if we want to reconstruct the mass dependence 
of fj,(uih)- We choose to approximate the functional form to be Gaussian, such that if a SM Higgs is present at mass 
m/j the number of signal events observed when applying the cuts, and hence searching, for a Higgs of mass m' h , is 

n s (m h ) = n s {m h )ex_p ^ ; (9) 

where a is the mass resolution for that search. For the case of multiple scalars the ATLAS h — > WW search has 
been studied in [2D]. The mass-dependence of the signal after cuts has been calculated and is shown in Fig 2. of [2DJ 
for hypothetical Higgs-like scalars of mass 125 and 170 GeV. One can see that due to the nature of the search the 
signal is not a delta-function centered at the Higgs mass but is rather a smooth distribution peaked at the true mass. 
We find that a Gaussian provides a good fit to the data, and so we assume that /i(m/i) follows the same functional 
dependence. 

Since the signal strength modifier is, by construction, normalized such that if a SM Higgs of mass mh is present in 
the data the strength modifier must be /j,(rnh) ~ 1, we normalize the Gaussian distribution to have a peak value of 
1. Given this assumption, combined with the approximate experimental resolution of the search channel, <7i(m), we 
estimate the strength modifier contributed by a SM Higgs of mass rrih in a particular search channel to be 

-tm—m ) 2 

H l (m) = exp 2-i 2 (4) , (10) 

where the normalization is chosen for a SM Higgs. Clearly, to extend this to a non-SM Higgs one includes depen- 
dence on any additional parameters by rescaling production cross-sections and branching ratios accordingly. Now to 
construct a likelihood for two Higgs scalars of mass and mi we estimate overlap of strength modifiers through the 
combination 

l i i {^m- h ) = ^) + ^)e , M >,m s ) = M |(u;) + /4(u,) e . (11) 

In this way, if the mass splittings far exceed the experimental resolution the strength modifiers become independent 
and the likelihood factorizes into individual likelihoods for the independent scalars. However as the masses approach 



The parameter dependence of the individual likelihoods is not necessarily independent. In the case considered here one might wish, for 
example, to increase the mixing, which increases the signal from s to explain some excess at mg 2> ra;. However doing so decreases the 
signal from h, which may be disfavored by the likelihood for h. 
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one another signal overlap becomes important, and in the limit where the masses are equal the strength modifiers 
simply add together, as expected. Alternatively one can think of this as the signal from one scalar acting as known 
background in the search for the other. This method is clearly approximate, however it should give a reasonable 
estimate of the combined likelihood given the available information and is useful to determine the mass range in 
which the factorized likelihoods can be trusted. 

We are combining multiple channels and so we must use different resolutions for each channel. CMS reports the 
approximate mass resolution of the individual channels in 41 , which we use, taking the maximum value whenever a 
range is quoted. We use 3% for the ZZ channel since this is the largest resolution in the individual ZZ sub-channels 
which are sensitive to a light Higgs. Our results will not be sensitive to this choice, since the dominant source of 
signal overlap is in the low-resolution channels unless mj k, m%. For ATLAS some resolutions are reported in |39j 
which are similar to those for the CMS searches. When not reported, we assume the same resolution as in the CMS 
searches. This assumption does introduce additional error whenever considering limits on scenarios where the scalars 
are close in mass, however in this region the dominant ATLAS sensitivity is in the h — > 77 and h — > ZZ channels, with 
published resolutions, and the h WW channel for which, by comparing with the results of [20], the assumption of 
20% mass resolution is valid. As a result it is likely that the overall error introduced into the combined limits and 
fits through this assumption is subdominant to other sources of error. The chosen mass resolutions are detailed in 
Table H 



TABLE I: Approximate Light Higgs Search Mass Resolutions 



Channel 


Resolution [%} 


77 


3 


ZZ 


3 


bb 


10 


TT 


20 


WW 


20 



C. Estimating the Importance of Signal Overlap 

To estimate the impact of signal overlap on best fit parameters we consider the case with pure Higgs-singlet mixing, 
setting a = 0. In this case strength modifiers for different search channels re-scale in the same way, simplifying the 
analysis. To perform this estimate we only use the 7 TeV data since best fit parameters and confidence contours are 
available for both CMS and ATLAS for the h — > bb, tt, WW, ZZ and h — » 77 search channels. 

In Fig. [T]we plot, in black, the best fit mixing angle as a function of the singlet mass, mg, for the combination of 



likelihoods of both scalars with signal overlap included according to Eq. (111. Due to a deficit in background events a 
negative Higgs event rate is, at some masses, preferred by the data. However, since we are fitting to a model restricted 
to real mixing angles, such negative event rates are not within the parameter space of the model, and these points 
usually correspond to a best-fit value of 9 = 0. We calculate 95% confidence bands by finding the mixing angle at 
which A(— 2 log A) = — 2(logA(/i) — logA(/i)) = 3.84. As argued in the caption, for < 210 GeV the overlap in 
signal clearly becomes important, and the simple product of likelihoods should not be used. 

However, since the best fit values shown in black do accommodate the signal overlap to some degree, we can still 
extract some qualitative features. Typically for masses mg < 200 GeV the preferred mixing angle is 9 « 0. This is 
due to two dominant effects. First of all, small mixing angles are preferred for the fit of = 125 GeV to the 7 
TeV data since the signal at this mass prefers /1 ~ 0.8, and larger mixing angles which reduce the signal further are 
penalized. Second, the strong limits for mj < 200 GeV also prefer the signal from the friend to be small, requiring a 
small mixing angle. One can also see that for mj the mixing angle essentially becomes unconstrained. This is 

due to the fact that in this case each strength modifier is almost independent of 9 since fi « sin 2 (#) + cos 2 (9) w 1. 



IV. RESULTS 



A. Data 



As demonstrated in Sec. IIIC for nig > 210 GeV we can effectively treat the likelihoods for the scalars individually, 
taking the product to find the combined likelihood. In this mass range the strength modifiers for the friend in the 



7 



1.0 
0.8 
0.6 
.9 0.4 
0.2 
0.0 





2011 Data 




a = 


if 















120 140 160 180 200 220 240 

m 

FIG. 1: The best fit mixing angle as a function of the Higgs friend mass, mg, for the combined likelihood with signal overlap 
included (black) and omitted (red). 95% confidence bands are also shown. Above ms ~ 210 GeV the difference between 
both methods becomes negligible, demonstrating that above this mass the simple product of individual likelihoods can be 
trusted. Below this mass the overlap of signal becomes important, suggesting that the simple individual likelihood products 
lose accuracy. 



relevant channels scale in the same way, since the diphoton search is not sensitive, and we only require the likelihood 
for the combined channels. Preliminary 7 and 8 TeV combined best fit values and confidence contours for the ATLAS 
searches were presented in [2] and so we employ these to construct the two-sided Gaussian likelihood. Preliminary 
7 and 8 TeV combined expected and observed 95% confidence limits for this mass range were presented for CMS in 
PQ, and so we use this data to estimate the combined likelihood in this region. In [37] it is shown that the best fit 
strength modifier can be simply approximated by the difference of the observed and expected upper limits, and the 
error is given by a « 

For the Higgs likelihood at 125 GeV we require information on individual channels since the parameter dependence 
differs for the 77 channel from the other channels. For this low mass region in [3j preliminary best fit values and 
confidence contours for the 77 and ZZ channels for combined 7 and 8 TeV runs at ATLAS are given. We estimate 
the full ATLAS likelihood function using these channels combined with 7 TeV values for the WW, bb and rr channels 
[59"] . all at 125 GeV. For CMS preliminary best fit values and uncertainties for the combined 7 and 8 TeV runs for all 
channels have been presented for a Higgs at 125 GeV pQ, and we use this to estimate the CMS likelihood function. 
Since these fits have been presented for a 125 GeV Higgs we assume this Higgs mass throughout. We also use the 
recent Tevatron data [22], taking the best fit values and uncertainties for the 77, WW and bb channels. 

With this information we can estimate the full likelihood for both scalars to determine whether a social Higgs allows 
for any improvement in fitting the data when compared with a SM (antisocial) Higgs. 



B. Higgs Friend Scenario 

First we consider the friend scenario, where the Higgs is mixed with a singlet scalar. On the left-hand panel of 
Fig. [2] we plot the best fit mixing angle as a function of the singlet mass for mg > 210 GeV. The best fit for the limit of 
vanishing invisible width is shown in black. We also plot results allowing for the decays BR(s — > 2h) following Eq. ^ 
with k = 0.5 in blue. Constraints are weakened by this effectively invisible width due to suppression of the signal, 
and fits previously requiring some mixing now require greater mixing due to dilution of the signal at high masses. 

On the whole, since the best fit points satisfy s'm 2 (8) < 0.15 and are consistent at 95% with sin 2 (6') = it is clear 
that the SM provides almost as good a fit as the mixed model, and around 400 GeV the Higgs friend scenario provides 
no advantage over the SM, even though some mixing is preferred for the 125 GeV signal. This is not true, however, 
at one point near ~ 210 GeV, where, due to a small excess in ZZ events, the Higgs friend scenario actually 
provides some improvement in fit over the SM, and the SM (0 = 0) actually lies on the 95% confidence contour. This 
improvement in fit is not great enough to suggest strong evidence for the presence of a Higgs friend, but is interesting 
nonetheless. If considered for just the singlet, this excess prefers \i ~ 0.15, and the Higgs signal at 125 GeV prefers 
fi ~ 0.92 and so both excesses fit the Higgs friend model well for sin 2 (8) ~ 0.15. Small excesses or weak limits at 



other masses also allow for a best fit with non-zero 8, although 8 = lies within the 95% confidence contour. 

Also, it is interesting to note that in some cases mixing angles as small as sin 2 8 ~ 0.1 are disfavored at the 95% 
level, showing that the Higgs searches are in some cases already sensitive to relatively minor modifications of the 
Higgs sector. We can also consider the case where the friend is effectively decoupled, with mass beyond the current 
sensitivity of the LHC searches. In this case we find that sin 2 (#) = 0.08 ± 0.19 for = 125 GeV, consistent with 
the SM. 
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FIG. 2: The best fit mixing angle as a function of the scalar mass, mg, at high masses. The limit of vanishing s — > 2h branching 
ratio is shown in black and for k = 0.5 (see Eq. |3|) in blue. 95% confidence contours are also shown (dashed) in corresponding 
colors. On the left panel we set a = and allow only mixing with the Higgs friend, with 95% confidence contours determined 
via A(— 2 log A) = 3.84. On the right panel we allow for enhanced h — ¥ 77 decays in the Higgs accomplice scenario and find 
the best fit values of 9 and a. In this case we find 95% confidence contours by finding the maximum value of 8 for which 
A(— 2 log A) = 5.99. On both panels it is clear that due to strong limits in the region 380 GeV < mg < 450 GeV the SM is 
preferred over both scenarios. However, due to mild excesses or weak limits at other masses both scenarios can slightly improve 
the fit to the data in comparison to the SM. Due to a small excess in h — > ZZ events in ATLAS and CMS, around rn% w 210 
GeV the accomplice scenario becomes marginally preferred over the SM at 95%, while for the friend scenario the SM lies on 
the 95% confidence contour. Following the discussion of Sec. |III C| the reader should keep in mind that in the region below 
mi ~ 210 GeV some error is introduced by neglecting signal overlap. 



C. Higgs Accomplice Scenario 

Whenever we allow for enhancement of the decays h — > 77 by coupling the friend to photons then, regardless of the 
value of 8, we can always choose the coupling, a, such that > 1 can be reproduced for the Higgs signals at 125 
GeV. The other search channels only constrain 8. Furthermore, as the diphoton searches look for resonances below 
150 GeV, for m s - > 150 GeV the likelihood function for s is independent of a. Hence a allows the freedom to achieve 
the desired Higgs diphoton rates without degrading the fit to the Higgs friend. For the 125 GeV excess there is a 
tantalizing hint that the 77 channel might be enhanced, while other rates might be suppressed, perhaps suggesting 
non-zero mixing angles. As such the Higgs accomplice scenario allows for a better fit to the data than the simple 
Higgs friend scenario. 

On the right-hand panel of Fig. [2] we show the best fit values for 8 as a function of mg. The 95% confidence contours 
are found by finding the largest possible value of 8 such that A(— 2 log A) = 5.99. Thus a best fit value of a, which is 
typically an 0(1) number (when positive), is concealed within this plot. 

For most of the mass range sin 2 (#) < 0.2 is preferred, with a slight increase in the best fit values due to the ability to 
accommodate the slightly suppressed rates in the non-diphoton channels at 125 GeV without suppressing the diphoton 
rates. At high masses weaker limits allow for larger signal, and hence mixing angles, to be accommodated for the 
Higgs accomplice while simultaneously fitting the slightly suppressed channels for the Higgs. The best fit value near 
210 GeV has also increased slightly in comparison with the friend scenario, as expected, and a slight bump in the 
lower 95% confidence limit shows that, at 95%, the SM is disfavored in comparison to the Higgs accomplice scenario 
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whenever rtiz = 125 GeV and wis — 210 GeV. Thus the Higgs accomplice scenario accommodates a better fit due to 
the additional source of 77 decays, however, the preference for this scenario is still not particularly strong. 
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FIG. 3: Best fit points and 68% and 95% confidence contours, corresponding to A(— 2 log A) = 2.28, 5.99 for the specific scenario 
of a Higgs accomplice beyond collider reach (left panel) and a Higgs accomplice at 210 GeV (right panel). The left panel shows 
the fit for the Higgs to solely the 125 GeV data since the accomplice is decoupled. Due to enhancement of h — > 77 and 
suppression of other channels, rather large mixing angles are preferred, alongside 0(1) values of a. The SM, sin 2 (0) = 0, is 
within the 95% confidence contour. The change in —2 log A, which can be interpreted as the change in x 2 , from the best fit 
point to the SM is also shown, and is greater than the number of extra parameters introduced. When signal from the Higgs 
accomplice is included at 210 GeV, and the likelihoods for both scalars are combined (right panel), the overall fit is improved 
significantly and the SM becomes marginally disfavored at greater than 95%. 

In Fig.[3]we show the best fit values as well as 68% and 95% confidence contours, corresponding to A(— 2 log A) = 2.28 
and A(— 2 log A) = 5.99, for a and 9 whenever s is decoupled and doesn't contribute any signal in the search window 
(left panel) and whenever signal from the friend is present at mj = 210 GeV (right panel). The left panel shows that 
with the friend decoupled the best fit points prefer substantial mixing, since the 77 rates can be fit independently of 
9. The SM at 9 = is within the 95% confidence contour, showing there is no strong preference for this scenario. 

Comparing both panels of Fig. |3j we see that the suppression of the non-diphoton event rates at 125 GeV requires 
larger mixing angles than are required to fit the small excess at 210 GeV, and the best fit mixing angle moves to 
smaller values when the friend is included. On the right-hand panel the SM, 9 — 0, lies outside the 95% confidence 
contour, demonstrating that if the accomplice mass is close to 210 GeV the Higgs accomplice scenario gives a definite 
improvement over the SM in fitting the data. In particular, a change of A(— 2 log A) ~ 8 between the SM and the 
best fit point indicates that the improvement in fit is not negligible, since only three new parameters are introduced 
(if one considers the accomplice mass to be fixed). 9 

Overall, the Higgs accomplice scenario does allow for a significant improvement in fitting the data over the SM, 
however this improvement is not great enough to claim support, or evidence, for a social Higgs. 

V. ELECTROWEAK PRECISION CONSTRAINTS 

As pointed out in [43 , if the Higgs mixes with a friend then precision electroweak observables are altered in 
comparison to the SM. In particular, W and Z boson couplings to the Higgs are suppressed, and the friend can also 



9 One might worry that the model is effectively over-fitting the data. Combining the measurements for the 125 GeV and 210 GeV Higgs 
searches one finds a log-likelihood value (which can be considered as approximately representing the \ 2 value) of —2 log A = 17 for the 
SM with D = 13 + 2 — 2 = 13 degrees of freedom (125 GeV channels, 210 GeV combined channels, and m^, /i). For the Higgs accomplice 
model, at the best fit parameters shown on the right-hand panel of Fig. [3] we find —2 log A = 9 for D = 13 + 2 — 2 — 3 = 10 degrees 
of freedom. Since X 2 / D ~ 1 for this model the data is not being over-fitted, and much of the improvement compared to the SM comes 
from fitting the excess at 210 GeV, which can not be accommodated within the SM. 
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enter at one loop into self-energy graphs. Here we study the differences in the S and T parameters [HI 25] relative 
to the SM for the Higgs friend model. We calculate these differences at one loop by taking the Higgs contributions 
to S and T from [151 SZ] and re-scaling them by cos 2 (9). We also add a similar contribution for the friend and then 
subtract off values for a SM Higgs at 125 GeV. 
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FIG. 4: Contours of AS = S(h,s,6) - S(h) and AT = T(h,s,6) - T(h), for the simple Higgs friend model. We have set 
= mi = 125 GeV. For the majority of parameter space this model is consistent with electroweak precision data at la. 



In Fig. [4] we show contours of the change in the S and T parameters relative to the SM with a Higgs at 125 GeV. 
In [48] for a Higgs mass in the range 115.5 < mu < 127 GeV the S and T parameters are given as S = O.OO^q \q 
and T = 0.02tg;^, so in this model deviations from the SM are typically within la, even with relatively large mixing 
angles, hence electroweak precision places no strong constraints on models with a relatively light Higgs friend. Heavier 
Higgs friends are less consistent with electroweak precision constraints, however agreement at close to la can still be 
found for friends with masses greater than 1 TeV |43j . 

Statements about precision electroweak observables are model-dependent and if additional electroweak-charged 
fields are present, as in the Higgs accomplice scenario, or as would be expected in a complete model which addresses 
fine-tuning issues, then further alterations to the S and T parameters would arise. In either case one must then consult 
the particular model to establish consistency with electroweak precision data. As such, the bounds shown here should 
be considered a demonstration of consistency in the friend scenario, rather than a reflection of the consistency of a 
possible underlying theory. 



VI. DISCUSSION 



The 125 GeV resonance recently discovered at the CERN LHC exhibits properties consistent with the SM Higgs 
boson. Only analyses of future data can convincingly determine whether or not it is indeed the SM Higgs, however, 
since there is currently no strong evidence to the contrary, it is now possible to constrain scenarios where the Higgs 
properties are significantly altered. Furthermore, null results in Higgs searches at other masses already place strong 
bounds on neutral scalars with Higgs-like production and decay properties. Motivated by this observation, in this 
work we have examined the impact this has on two simple models, the Higgs friend and Higgs accomplice scenarios, 
which we believe may act as simplified models for theoretically motivated extended Higgs sectors. Both scenarios 
are still compatible with the data, however large mixing angles 6 > 7r/4 are typically disfavored at the 95% level. 
Small mixing angles satisfying sin 2 ((9) < 0.2, can improve the overall fit for the Higgs at 125 GeV, especially if the 
model accommodates enhanced couplings to photons. However, the improvement in fit is, in the majority of cases, 
not statistically significant. The only exception is for a Higgs accomplice with a mass near 210 GeV which allows for 
an improvement in fitting the data at greater than 95% significance when compared to the SM . However, keeping the 
numerous uncertainties and the low statistical significance of the excess in mind, one cannot interpret this as evidence 
for a bone-fide Higgs friend or accomplice. 

Note added: While this work was being completed a similar analysis was completed [35], which considers the 
Higgs accomplice scenario whenever the friend is decoupled. In [35] CMS and ATLAS results for h —> tt searches 
and Tevatron results for h — > WW and h — > 77 are not included. This leads to smaller preferred mixing angles than 
we find here. Omitting these channels from our analysis we find results in agreement with those presented in [35) . 
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